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1. Introduction 

The concern in the mid 1970s regarding health 
effects from exposure to electric and magnetic 
fields in the vicinity of power lines has shifted in re- 
cent years to health effect concerns from exposure 
to power frequency magnetic fields in residences, 
the work place, and in transportation systems [1-3]. 
The magnetic fields in these environments can be 
highly nonuniform, particularly near electrical 
equipment such as motors, transformers and heat- 
ing elements. This paper considers the difference 
between the calculated average magnetic flux den- 
sity, 5av, as determined using magnetic field meters 
with single-axis and three-axis circular coil probes, 
and the calculated magnetic flux density at the cen- 
ter of the probes. Bo assuming the field is produced 
by a small loop of alternating current, i.e., a mag- 
netic dipole. The magnetic dipole field is chosen as 
the relevant field because to a good approximation 
its geometry simulates the field geometry of many 
electrical appliances and equipment [4]. The differ- 
ence between Bm and Bo can be regarded as a 
source of measurement uncertainty because the 



center of the probe is normally considered the mea- 
surement location. While differences between B^ 
and Bo will be small in many situations, e.g., near 
ground level in the vicinity of power lines where the 
field changes slowly, the difference can become sig- 
nificant in the highly nonuniform magnetic fields 
close to electrical equipment. 

In this paper, two comparisons are made: (1) the 
maximum average magnetic field determined using 
a single-axis probe, 5avi, with Bo as a function of 
r/a where r is the distance between the magnetic 
dipole and the center of the probe, and a is the ra- 
dius of the probe, and (2) the resultant magnetic 
field determined using a three-axis probe with Bo as 
a function of r/a . The resultant magnetic field, 5av3, 
is defined as [5] 



B^y3 = VBl+Bl+Bl, 



(1) 



where Bi, Bi, and Bz are average magnetic field 
components as measured by three orthogonally 
oriented coil probes. 
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Comparison (1) is made because maximum mag- 
netic field values are sometimes measured, using 
single-axis field meters, to characterize the mag- 
netic field [5,6]. However, for a given value of ria , 
it will be seen that the difference between 5avi and 
Ba will be a function of the orientation of the mag- 
netic dipole relative to the probe. Because the rela- 
tive orientation is not known during most 
measurements, what is of interest is the largest dif- 
ference between Bm\ and Bo for a given value of rIa . 
This largest difference will be designated ABmai. 

The quantity ABmaxi is determined in the follow- 
ing way. The single-axis probe is rotated for fixed 
values of rIa and the spherical coordinate, 6, 
[Fig. 1(a)] until the largest average magnetic field, 
Bw\, is found. This value of 5avi is compared with 
the magnetic field at the center of the probe, Bo 
and the difference is recorded. The orientation of 
the magnetic dipole with respect to the probe is 
then varied by moving the probe to another loca- 
tion while keeping rIa fixed, i.e., by changing 6 in 
Fig. 1(a). The probe is rotated again until the 
largest average magnetic field, Bmi, is found. Bwi is 
again compared with the magnetic field at the cen- 
ter of the probe. Bo, and the difference is recorded. 
This process is repeated for other dipole orienta- 
tions (i.e., angle 6) until the largest difference, 
ABmaxi, is found. An example of this process is 
shown in Sec. 3.1. 



single-axis 
probe 




ioop of 

aitemating 

current 



(a) 



(b) 



Fig. 1. (a) Single-axis and (b) three-axis circular coil probes in 
dipole magnetic field produced by small loop of current. 

Comparison (2) is made as a three-axis probe is 
rotated about three axes parallel to the three 
Cartesian coordinates x, y, and z. The difference 
between Bavs and Bo will vary as a function of rota- 
tion angle, but what will be of interest again is the 
largest difference, ABmaxs, for a given rIa. Also as 
for comparison (1), because the relative orienta- 
tions of the magnetic dipole and the three-axis 
probe will be unknown in most measurement situa- 
tions, Bav3 will be examined as a function of rIa and 



the spherical coordinate, 6, in order to determine 
the largest difference, ABmaxs. 

2. Expressions for Average Magnetic 
Flux Density 

In the derivations given below, it is assumed that 
the cross sectional area of the wire in the coil 
probes and the opposing magnetic field produced 
by current induced in the probes are negligible. In 
addition, we assume for the three-axis probe that 
the three orthogonally oriented coils have circular 
cross sections of equal area. These assumptions ei- 
ther can be met in practice or can be taken into 
account via a calibration process. 

2.1 Single-Axis Probe 

The average magnetic flux density, Bav, for a sin- 
gle coil probe with cross sectional area A is given 
by 



Bav=j^ J J B-hAA, 



(2) 



where 6A is an element of probe area, n is a unit 
vector perpendicular to A, and B is the magnetic 
flux density. In spherical coordinates, the magnetic 
flux density for a small current loop of radius b is [7] 



B = 



^ 3 cos Our + ^^3 sm 9 Ue , 



2r' 



4^3 



(3) 



where no is the permeability of vacuum, / is the al- 
ternating current, and Ur and Ue are unit vectors in 
the directions of increasing r and 6, respectively. 
The assumption is made that b < <r and the sinu- 
soidal time dependence has been suppressed. The 
magnitude of the vector B given by Eq. (3) is Bo. 

For our purposes, it is convenient to express B in 
terms of Cartesian coordinates. This is accom- 
plished by using the following relations between 
spherical and Cartesian unit vectors and coordi- 
nates [8] in Eq. (3): 

«r=?sin 0COS0 +/sin 0sin<^+^cos 9 
UB=i cos cos +/ cos sin —k sin 9 
x=r sin 9 cos 
y =r sin 9 sin 
z=r cos 9 . 



(4) 



After some algebra, B can be expressed as 



B=i 



.3Cxz 



Ir' 



.;^ ^^1,(^-1). (5) 
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where r = v'jc^+y^+z^ and C is the constant 
iUo/6'/2. 

To obtain an expression for J?av, we consider 
without loss of generality a probe with its center at 
x=xa,y=0, and z=za as shown in Fig. 2. We 
restrict the orientation of the probe so that its area 
is bisected by the x-z plane and first consider rota- 
tions of the probe about an axis parallel to the 
y-axiSi i.e., they '-axis shown in Fig. 2. For these 
conditions, the area of the coil probe. A, will be 
part of the surface given by the equation 



z=ma(x-JCo)+zo, 



(6) 



where a is the angle of rotation, nia = tan a, 
XQ=r sin 0, and zo=r cos 0. The rotation of the 
probe corresponds to the rotation of this surface 
about the y'-axis, i.e., changing the slope of the 
surface (wia) described by Eq. (6). It should be 
noted that the angle of rotation, a, shown in Fig. 2 
is in the negative direction. 




surface: 
m» =tanQr 



*-y 



x=Xg- acosa 
x=Xg + acosa 



Fig. 2. Circular coil probe shown as part of a surface described 
by the equation z=m„(x-Xo)+Zii. The rotation of the probe 
corresponds to changing the slope of the surface, m„. The pro- 
jection of the probe cross sectional area onto the x-y plane will 
be an ellipse for av^O. The range of a is -90°<a< 90°. 



The unit vector perpendicular to the probe sur- 
face, n , is found by first taking the gradient [9] of 
the surface given by Eq. (6), VF(x,z), where 
F(x,z)=z-ma(x-Xo)-Zo and normalizing it to 
unit value. This leads to 



/J = ( -mj +k)/y/mi+ 1. 



The element of area, dA , is [10] 



(7) 



dA 



= ^(|)'+(&)'+ldxdy = Vm„^ + l dxdy. 



(8) 



where dx dy is an element of area in the jc-y plane 
bounded by the projection of the probe cross sec- 
tion onto the x-y plane (Fig.2). 

Combining Eqs. (2), (5), (7), and (8), the expres- 
sion for Bm becomes 



Bm — 



lira- 



H{-^-^M¥ -')}■-<'- 



(9) 



By substituting forz in Eq. (9) using Eq. (6), the in- 
tegrand becomes a function of jc andy. The integra- 
tion is first carried out analytically [11] over the 
variable y with (from Fig.2) 



- y/a^ - ((jc - JCo)/cos of ^y. 



yja^ - ((x - JCb)/cos of 



The resulting expression for B^ is 



n --- ^ 



J dxipioXP ■VP^') 



Ja^-( x-Xo \2 
* \cos a/ 



(,.+e.)(,.+e.+,._(|i^)3)- 



'* Vcos a/ 



(,.+e.).(,.+e.+,._(-^).) 



:l- 






dx- 



Ja^-f x-Xo V 
^ Vcos a) 



(-'■'H^'-^Q'-'-'-i^yy' 



(10) 
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where P = (zo-/M,;co) and Q^=(mjc + py. 

The integration over x is then performed 
numerically using Simpson's Rule with the limits of 
integration given by (Fig. 2) 

Xo—a cosa^x^Xo+a cos a, 

where a is restricted to - 90° < a < 90°. 

Bm is evaluated for fixed values of 6 and r/a as a 
is varied until a maximum average flux density, 5avi, 
is found. 5avi is then compared with Bo. As noted 
earlier, the process is repeated for the same r/a but 
different values of 6 until the largest difference, 
ABrasau is determined. Because we are seeking the 
maximum value of B^, we do not consider further 
rotations of the probe because once B^i is found, 
additional rotations are expected to lead to smaller 
values of 5av. This is most readily seen at 6 equal to 
0° and 90° for all values of r/a . Bmi occurs at a = 0° 
and rotating the probe further results in smaller 
values of 5av. 

2.2 Three-Axis Probe 

In this section, expressions are developed for the 
average magnetic flux density for each coil of a 
three-axis probe as the probe is rotated about axes 
which are parallel to the x-, y-, and z-axes. After- 
wards, for fixed values of 6 and r/a, the average 
magnetic field values from the three orthogonally 
oriented probes are combined according to Eq. (1) 
to obtain 5av3 which is then compared with Bo. As 
before, the process is repeated for different values 
of 6 until the maximum difference between Bm3 and 
Bo, ABmai, is found. 

It is noted that combinations of rotations about 
the different axes will not be possible using the ex- 
pressions that are developed. That is, it will not be 
possible to calculate 5av3 following rotations about 
two or three axes. This represents a limitation on 
the results and prevents us from learning whether 
there are significant effects on the value of ABmaxi 
due to multiple rotations. Nevertheless, the depar- 
tures from Bo that are determined from rotations 
about each of the three axes will let us know what 
differences are possible as a function of r/a . 

Figure 3 shows the orientation of the probe with 
respect to the magnetic dipole before rotations 
about each axis are considered. As in the previous 
section, the center of the probe is located at x =xo, 
y=0, and z =Zo. The individual probes are labelled 
PI, P2, and P3, and initially PI lies in the x'-y' 
plane, P2 is in tht y'-z' plane, and P3 is in thex'-z' 
plane. When rotations are performed about the 
x'-axis (rotation angle jS), the angle that P2 makes 



with the magnetic field remains unchanged. Simi- 
larly, rotations about they '-axis (rotation angle a) 
and z '-axis (rotation angle y) leave P3 and PI, re- 
spectively, "unchanged." This means that for con- 
stant values of r/a and 6, the average flux density 
values for these "fixed" probes remain constant as 
rotations of the probe occur. 




Fig. 3. Geometry for three-axis probe with center of probe in 
x-z plane. Varying the angle a result in rotations of probes PI 
and P2 about the ^'-axis ("a rotations") while orientation of 
probe P3 with respect to the dipole remains unchanged. Varying 
the angle /3 results in rotations of probes PI and P3 about the 
a: '-axis ("/3 rotations") while orientation of probe P2 remains 
unchanged. Vaiying the angle y results in rotations of probes P2 
and P3 about the z'-axis ("y rotations") while orientation of 
probe PI remains unchanged. 

We begin the derivation for the three-axis probe 
by noting that part of the problem has already been 
solved in the previous section. That is, the expres- 
sion for 5av following rotations about the y'-axis 
(a rotations) is given by Eq. (10). This expression is 
used to calculate the average flux density for 
probes PI and P2 by considering pairs of the angle 
a which differ by 90°. The average flux density for 
the third probe, P3, is zero for this case because no 
component of the magnetic field is perpendicular 
to the area of the probe for any value of a or 6. 

The derivations for B^ following /3 or y rotations 
parallel the derivation for the a rotations. Examin- 
ing the case of j8 rotations first, and referring to 
Fig. 4, it can be seen that the probe area, A , is part 
of the surface given by the equation 



z=m^>',+Zo, 



(11) 
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surface: 







-T35 



Fig. 4. Circular coil probe stiown as part of a surface described 
by the equation z=mpy+zo. The rotation of the probe corre- 
sponds to changing the slope of the surface, ntp. The projection 
of the probe cross sectional area onto the x-y plane will be an 
ellipse for j8;*0. The range of iS is -90°< /S < 90°. 



where m^ = tan ^ and zo=r cos 9. The rotation of 
the probe corresponds to rotation of this surface 
about the jc'-axis or alternatively, changing its 
slope, ntp. The unit vector normal to this surface is 



/i=(-m^+A:)/\/m|+l, 
and the element of area, dui , is 



(12) 



dA 



■■y/[^J+[f)l\l<ixdy=Vmj + l dxdy. 



(13) 



Combining Eqs. (2), (5), (12), and (13), the expres- 
sion for the average flux density following )3 rota- 
tions is 



B^p=-^ I \ B-ndA = 



Jy Jx 



3Cm, 
' Ittu' 



\.\3^'y 



lira 



[[[^-7)^^y- (14) 



Substituting for z in Eq. (14) using Eq. (11), the 
integrands become, recalling that r = y/x^+y^+z\ 
a function of jc andy. The integration is first car- 
ried out analytically over the variable x with (Fig. 4) 



Xo-Va^-iy/cos pf^x^xo + Va^-(y /cos py. 

The y -integration is performed numerically with 
the limits of integration given by (Fig. 4) 

—a cos p^y^ a cos )3, 

where )3 is restricted to -90° < /3 < 90°. 

Equation (14) is used to calculate the average 
flux density for probes PI and P3 by considering 
pairs of the angle p which differ by 90°. The aver- 
age flux density from the remaining probe, P2, re- 
mains constant during the )3 rotations and is 
determined from the expression for average flux 
density following y rotation (■y=0)— which is de- 
veloped below. The reader is cautioned that during 
the numerical integration over the variable y, the 
denominator in the integrand vanishes for y=0 
when = 90°. 

In deriving the expression for average flux den- 
sity following y rotations, we note as shown in 
Fig. 5 that the probe area, A , is part of the surface 
given by the equation 



x=myy +xo. 



(15) 



where my = tan-y and xo=r sin 9. The rotation of 
the probe corresponds to rotation of this surface 
about the z'-axis, i.e., changing the slope of the 
surface, my. The unit vector normal to this surface 
is 



n=(i-m y])l\/my + 1, 



(16) 
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and the element of area, cL4 , is 



-W^' 



+ 1 d>' dz = ■\/my + 1 d>' dz. 
(17) 



TOP VIEW 



surface: 
my = tan 7 




Fig. 5. Circular coil probe shown as part of a surface described 
by the equation x =myy +xo. The rotation of the probe corre- 
sponds to changing the slope of the surface, my. The projection 
of the probe cross sectional area onto the y-z plane will be an 
ellipse for ypsO. The range of y is -90° <>' < 90°. 

It should be noted that the projection of the probe's 
cross sectional area, unlike the previous two cases, 
is onto the>'-z plane for y rotations (the projections 
for a and /3 rotations were onto theA;->' plane) and 
this fact affects the partial derivatives in the expres- 
sion for dA , Eq. (17). 

From Eqs. (2), (5), (16), and (17), the expression 
for average flux density following y rotations is 



^-^=^11 ^■"'^^ 



JC_ 

lira- 



IKP- 



^)dzd. 



(18) 



By substituting for a: in Eq. (18) using Eq. (15), the 
integrand becomes a function of >> and z. The inte- 
gration is first carried out analytically over the vari- 
able z with (Fig. 5) 



Zo-\/a^-(y/cos-y)^ ^z ^Zo+\/fl^-(y/cos-y)l 

They -integration is performed numerically with the 
limits of integration given by 

- fl cos-y ^>' ^ fl cos-y, 

where y is restricted to -90°<'y<90°. 

Equation (18) is used to calculate the average 
flux density for probes P2 and P3 by considering 
pairs of the angle, y, which differ by 90°. The aver- 
age flux density from the remaining probe, PI, 
remains constant during the y rotations and is 
calculated from the expression for average flux den- 
sity for a rotations [Eq. (10)] with a set equal to 
zero. 

3. Results of Calculations 
3.1 ilfinuxi for Single-Axis Probe 

Using Eq. (10) and following the procedure de- 
scribed after Eq. (1), values of the maximum aver- 
age magnetic field, B^i for fixed values of r/a and 9 
were calculated and compared with the correspond- 
ing value of Bo. Figure 6 shows the differences in 
percent between Bmi and Bo for r/a = 3 and for rep- 
resentative values of 9 between 0° and 90° (because 
of symmetry arguments, one can infer the corre- 
sponding percentages for 9 between 90° and 180°). 
The largest difference, ABmaxi, is -14.6% and 
occurs when the single-axis probe is located along 
the axis of the magnetic dipole, i.e., the z-axis. The 
negative difference between Bavi and Bo decreases 
as 9 increases and turns positive near 9 = 90°. This 
pattern also occurs for other values of r/a greater 
than 3. Figure 6 also shows the largest negative and 
positive differences in percent for r/a equal to 5, 8, 
10, and 12. The largest negative differences must be 
considered part of the measurement uncertainty 
when the probe-dipole geometry is unknown, which 
will be the case for example when magnetic field 
measurements are performed near many appli- 
ances. A tabulation of ABmau as a function of r/a is 
given in Table 1. 

The calculations are not carried out for large val- 
ues of r/a because the accura(^ requirements for 
magnetic field measurements near appliances and 
other electrical equipment either have not been set 
or are not great. For example, the uncertainty 
tentatively allowed during calibration of magnetic 
field meters used for measuring magnetic fields 
near visual display terminals is ±5% [12]. 
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e = i5' 



e = 30° 



e = 45° 



e = 6o° 




e = 75° 



> e = 90° 



Fig. 6. Differences between values of B^i and Bo, in percent, 
for different locations of single-axis probe relative to magnetic 
dipole which is aligned along the z-axis. For a given value of r/a , 
where r/a is 5 3, the largest difference, ABmaxi, is negative and 
occurs when the probe is located along the z-axis. 



Table 1. Values of 4Bmaxi (single-axis probe) and 4Bmax3 
(three-axis probe) as a function of normalized distance (r/a) 
from magnetic dipole 



r/a 



AB„ 



i(%) 



4B™ari(%) 



3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 



-14.6 
-8.7 
-5.7 
-4.0 
-3.0 
-2.3 
-1.8 
-1.5 
-1.2 
-1.0 
-0.9 
-0.8 
-0.7 



-19.6 
-10.8 
-6.9 
-4.8 
-3.5 
-2.7 
-2.1 
-1.7 
-1.4 
-1.2 
-1.0 
-0.9 
-0.8 



3.2 J0m»c3 for Three-Axis Probe 

The differences between 5av3 and Bo are consid- 
ered in three steps. First, values of 8^3 are calcu- 
lated following a rotations using Eq. (10) and 
compared with Bo for fixed values of r/a and repre- 
sentative values of 6 between 0° and 90°. As noted 
earlier, the largest difference between 5av3 and Bo at 
each point, a "local maximum difference," is 
recorded. In the discussion that follows, the "local 
maximum difference," will be referred to simply as 



the "difference." The above procedure is repeated 
for P and y rotations. 

The differences between /3av3 and )3b following a 
rotations are plotted in Fig. 7 for r/a equal to 3, 5, 
8, 10, and 12. The numbers in Fig. 7 represent dif- 
ferences in percent. The pattern observed for all 
values of r/a is that the difference at a point follow- 
ing a rotation is always negative and becomes more 
negative as 6 increases to 90°. 



e = i5' 



e = 30° 



e = 60° 




e = 75° 



> e = 90° 



Fig. 7. Differences between values of Bav3 and Bo in percent, for 
different locations of three-axis probe relative to magnetic 
dipole, following a rotations. The differences are always negative 
and the greatest difference following a rotations occurs for 
= 90°. 

When the difference calculations are performed 
for 13 rotations using Eq. (14) and Eq. (18) (y set 
equal to zero), a different pattern emerges. The 
differences between Bm3 and Bo are observed to 
change in sign at different points as shown in 
Fig. 8. By comparing the results in Figs. 7 and 8, it 
can be seen that, except for 6 = 0, the differences 
following j8 rotations are all less than the corre- 
sponding (i.e., same r/a and 6 values) differences 
following a rotations. Although not indicated in 
Fig. 8 calculations show that for a given value of r/a 
>3, the differences between Bm3 and Bo following /3 
rotations for all values of 6 are less than the differ- 
ence following a rotations when 6 = 90°. 

Similar results occur following y rotations when 
5av3 is calculated using Eq. (18) and Eq. (10) (a set 
equal to zero) and is compared with Bo. Once again 
the differences between 5av3 and Bo change in sign 
as shown in Fig. 9. Also, by comparing results in 
Figs. 7 and 9, it is seen that, except for 6 = 0°, the 
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8 = 0' 



8 = 30° 



8 = 60° 




8 = 75° 



> 8 = 89.9° 



Fig. 8. Differences between values of B„3 and Bo in percent, 
for different locations of three-axis probe relative to magnetic 
dipole, following p rotations. The differences vary in sign de- 
pending on angle 0. For = 15°, there are several cases (indi- 
cated with the ± sign) for which the largest negative and 
positive differences are equal in magnitude. 



8 = 0° 



8 = 30° 



8 = 60° 




8 = 75° 



> 8 = 90° 



Fig. 9. Differences between values of S,v3 and So in percent, 
for different locations of three-axis probe relative to magnetic 
dipole, following y rotations. The differences vary in sign 
depending on angle 0. 

differences following y rotations are all less than 
the corresponding differences following a rotations. 
As for the case of /3 rotations, calculations show 
that for a given value of r/a ^ 3, the largest differ- 
ence following y rotations will always be less than 



the difference following a rotations when 6 = 90°. 
Therefore, of the three types of three-axis probe 
rotations considered, the greatest difference 
between 8^3 and Bq can be found following a rota- 
tions and ABaaia occurs when 0=90° for a given 
value of r/a . Table 1 provides a listing of ABnaoa 
values as a function of r/a . 

4. Discussion of Results 

Once it has been decided what constitutes an ac- 
ceptable level of uncertainty during magnetic field 
measurements near electrical equipment, the infor- 
mation in Table 1 should be considered when taking 
into account the various sources of measurement 
uncertainty. For example, if maximum magnetic 
fields at a distance r from appliances are to be mea- 
sured with a total uncertainty of less than ± 10%, 
magnetic Held meters with probes having radii a 
such that r/a ~3 would immediately be considered 
unsuitable. Field meters with single-axis probes 
having radii such that r/a = 5 would be suitable if all 
other sources of uncertainty (e.g., calibration pro- 
cess, freque ncy respo nse) amounted to about 8% or 
less, i.e., V5.7^+8^ = 9.8, where 5.7 is taken from 
Table 1 for r/a =5. 

The measurement uncertainties associated with 
using three-axis probes are less clear because we 
have considered only separate rotations about three 
axes to obtain the values of /IBmaxs . The percentage 
differences in Table 1 indicate what uncertainties 
can occur but they may not be the largest uncertain- 
ties due to the averaging effects of the probe. How- 
ever, until calculations can be devised which 
consider more complex rotations of three-axis 
probes, the ABmaxi values in Table 1 can serve as a 
rough guide when deciding what are acceptable 
probe dimensions. 
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